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ABSTRACT 

The present-day formation of cataclysmic variables (CVs) with brown dwarf 
(BD) secondaries (0.013 < < 0.075 M©) is investigated using a popula- 
tion synthesis technique. Results from the latest, detailed models for BDs have 
been incorporated into the population synthesis code. The present-day orbital 
period distribution of zero-age CVs (ZACVs) that form with BD secondaries is 
calculated. For our models, we find that ZACVs with BD secondaries have orbital 
periods in the range 46 min to 2.5 hrs. We also find that ZACVs with BD secon- 
daries comprise 18% of the total, present-day ZACV population. In addition, we 
find that 80% of ZACVs with BD secondaries have orbital periods less than 78 
minutes. This implies that 15% of the present-day ZACV population should have 
orbital periods shorter than the observed orbital period minimum for CVs. We 
also investigate the dependence of the present-day formation rate of CVs with 
BD secondaries on the assumed value of the common envelope efficiency param- 
eter, acEi for three different assumed mass ratio distributions in ZAMS binaries. 
Surprisingly, we find that the common envelope process must be extremely in- 
efficient {acE < 0.1) in order for CVs with BD secondaries not to be formed. 
Finally, we find that the progenitor binaries of ZACVs with BD secondaries have 
ZAMS orbital separations < 3 AU and ZAMS primary masses between ~ 1-10 
Mq, with ~ 75% of the primary masses less than ~ 1.6 Mq. Interestingly, these 
ranges in orbital separation and primary mass place the majority of the progeni- 
tor binaries within the so-called "brown dwarf desert. " The implications of these 
results are discussed in the context of V485 Cen and IRXS J232953.9+062814, 
the only known CVs with orbital periods shorter than the period minimum, CVs 
above the period minimum that are strongly suspected of containing BD secon- 
daries (such as WZ Sge), common envelope evolution involving very low mass 
secondaries, and the BD mass function in ZAMS binaries. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: evolution — 
stars: formation — stars: fundamental parameters — stars: low-mass, brown dwarfs — 
stars: luminosity function, mass function 
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1. Introduction 

The formation of cataclysmic variables (CVs) has been studied in detail by several 
authors (Politano 1988, 1996; deKool 1992; Han et al. 1995; Howell et al. 2001). However, 
there has been only one detailed study of the formation of CVs with brown dwarf (BD) 
secondaries to date, that of Politano (1988, 1996). This study used (by now) outdated physics 
for the radius-mass relationship of the BD and assumed that the primordial distributions of 
main-sequence (MS) -|- BD binaries was the same as the primordial distribution of MS -|- MS 
binaries. A more recent study considered the secular evolution of CVs with BD secondaries 
using up-to-date models for the brown dwarfs (Kolb & Baraffe 1999). However, this study 
did not calculate the formation of CVs with brown dwarf secondaries, but rather calculated 
the CV birthrate as a function of secondary mass for MS stars and extrapolated it into the 
BD regime. 

Significant advances both in our theoretical and observational understanding of brown 
dwarfs have occurred in the past five years (e.g.. Burrows et al. 1997; Chabrier & Baraffe 
1997, 2000; Kirkpatrick et al. 2000; Zapatero Osorio et al. 2000; Bcjar et al. 2001; Burgasser 
et al. 2002; Dahn et al. 2002). Such important advances now offer an excellent opportunity 
to apply these models to the low mass secondaries in CVs. We have incorporated into our 
population synthesis code (Politano 1996) current models of low-mass stars and BDs from 
the Lyon group. Their evolutionary code includes improved internal physics (Chabrier & 
Baraffe 1997, 2000) and outer boundary conditions based on non-gray atmosphere models 
(see Allard et al. 1997), appropriate to describe the complex structure of low mass stars and 
brown dwarfs. These improvements in the present theory have led to a better understanding 
of the properties of M-dwarfs of solar metallicity (Baraffe et al. 1998) and yield excellent 
agreement with globular cluster lower main sequences (Baraffe et al. 1997). The models 
extended to the brown dwarf regime (Chabrier et al. 2000; Baraffe et al. 2003) successfully 
reproduce the near-IR photometric and spectroscopic properties of the cool brown dwarf GL 
229B (Allard et al. 1997), as well as the young and hot brown dwarfs observed recently in 
young clusters (cf. Zapatero Osorio et al. 1997, 2000; Martin et al. 1998; Bejar et al. 2001). 

Further motivation for the present work is provided by the existence of two unusual 
systems: V485 Cen, a CV with an orbital period of 59 min (Augusteijn et al. 1993, 1996; 
Olech 1997), and IRXS J232953.9+062814, a CV with an orbital period of 64 min (hereafter 
RX2329+06; Jingyao et al. 1998; Uemura et al. 2001, 2002; Thorstensen et al. 2002). In the 
standard model of CV secular evolution, a CV evolves to shorter orbital periods, reaches a 
minimum period, and then begins to evolve back toward longer orbital periods (Paczyhski 
& Sienkiewicz 1981; Rappaport, Joss & Webbink 1982; see also King 1988 for a review). 
The orbital periods in both V485 Cen and RX2329-I-06 are shorter than either the observed 



-3- 



(78 min) or theoretically-predicted (~ 70 min) orbital period minimum (cf. Kolb & Baraffe 
1999) for CVs that contain relatively uncvolved main sequence secondary stars, and thus 
cannot be explained within the standard model of CV secular evolution. 

There are objects related to CVs with orbital periods less than 78 min, namely the AM 
CVn stars (Porb ~ 20-50 min; Ritter & Kolb 2003; Nelemans et al. 2001). AM CVn stars are 
believed to contain a white dwarf and a degenerate helium white dwarf or a semi- degenerate 
helium star (see Ulla 1994 for a review; also Nelemans et al. 2001 and references therein). 
This model is consistent with the notable very weak hydrogen lines in the spectra of these 
stars. However, both V485 Ccn and RX2329+06 have orbital periods longer than known 
AM CVn stars, and both systems contain hydrogen lines in their spectra (Augusteijn et al. 
1993, 1996; Thorstensen et al. 2002). The presence of hydrogen hues in the spectra of V485 
Cen and RX2329-I-06 thus precludes a model for these systems in which the secondary is a 
He star. 

It has been suggested instead that both V485 Cen and RX2329-I-06 contain a secondary 
in which a significant fraction of hydrogen in the core had been processed into helium prior to 
mass transfer (Augusteijn et al. 1996; Schenker & King 2002; Thorstensen et al. 2002; Uemura 
et al. 2002). By the time such a system reaches short orbital periods, significant erosion of 
the secondary will have exposed the He-rich core. Calculations of the secular evolution of 
such CVs predict that these systems can reach orbital periods as short as those observed 
in V485 Cen and RX2329+06 (Schenker & King 2002; Thorstensen et al. 2002). However, 
Politano (1988, 1996) also found that CVs that form with substellar secondaries could have 
orbital periods, at birth, as short as 34 min. Since his calculations relied on outdated 
physics for the substellar secondaries, his result needs to be confirmed using population 
synthesis calculations that contain up-to date models for the substellar secondaries in order 
to determine whether a CV that formed with a BD secondary is also a possible model for 
V485 Cen and RX2329+06. 

We also examine in this paper the dependence of the formation rate of CVs with BD 
secondaries on the assumed value of the common envelope efficiency parameter, acE- In the 
absence of readily incorporable results from detailed hydrodynamical models of the common 
envelope (CE) phase (cf. Iben & Livio 1993; Taam & Sandquist 2000 for recent reviews), 
population synthesis calculations have relied on a simple energy argument to relate the 
prc-CE and post-CE orbital parameters. Essentially, this argument states that when some 
fraction of the orbital energy released in the spiral-in process equals the initial binding energy 
of the primary's envelope, the CE becomes unbound and is presumably (somehow) ejected 
(Tutukov & Yungleson 1979). Although the exact prescription varies from author to author, 
a typical expression used in population synthesis calculations is given below (cf. Iben & 



-4- 



Tutukov 1985; di Stefano & Rappaport 1994; Politano 1996): 



OiCE{ 




JlGMeMe_j2GMeMc 

Ri Ri 



(1) 



where Mp is the total primary mass at the onset of the CE phase, Mg is the mass of the 
secondary, Mp is the core mass of the primary. Me is the envelope mass of the primary (Me 
= Mp-Mc), Ri is the radius of the primary at the onset of the CE phase (and is therefore 
equal to the radius of the primary's Roche lobe), ai and af are the pre-and post-CE orbital 
separations, respectively, 71 and 72 are dimensionless factors related to the structure of the 
giant star, and ac^; is a parameter that measures the efficiency at which orbital energy is 
able to be used in unbinding the CE. 

Prom equation 1, we see that the orbital energy depends linearly on the mass of the sec- 
ondary. Because of the small secondary mass, in considering the spiral-in of a BD secondary 
into the giant envelope, there may simply not be enough of an energy reservoir available in 
the orbit to unbind the envelope if the transfer of energy from the orbit is too inefficient. 
Instead the CE phase may result in merger of the component stars. While Politano (1996) 
found that CVs with BD secondaries were able to be formed, he assumed that the transfer 
of energy was 100% efficient (i.e., acE — 1)- This has been a fairly standard assumption 
in population synthesis studies of CVs (e.g., deKool 1992; Yungelson et al. 1995; Han et 
al. 1995; Howell et al. 2001). However, some detailed hydrodynamical calculations of the 
common envelope phase suggest that the value of acE is lower (~ 0.3-0.6; e.g., Iben & Livio 
1993; Taam & Sandquist 2000). We should note that detailed hydrodynamical calculations 
cannot follow the CE phase completely through envelope ejection as yet, so any estimates of 
QcE arc based on extrapolations from these calculations, and should be regarded as sugges- 
tive rather than definitive. The efficiency of the transfer of orbital energy into the common 
envelope remains a very uncertain quantity. 

In this paper, therefore, we quantitatively study the dependence of the formation rate of 
CVs with brown dwarf secondaries on the value of ace- In doing so, a question of particular 
interest is whether there is a minimum acE below which CVs with BD secondaries are not 
able to be formed. We note that Kolb (1993) calculated model populations of secularly- 
evolved CVs using two different values of the efficiency parameter, acE — 1 and acE — 0.3, 
and found little difference in the model populations. However, his study did not consider 
CVs that contained substellar secondaries at birth. In addition, his study did not consider 
values of acE less than 0.3. 

To the best of our knowledge, the only hydrodynamical calculations of the spiral-in 
of a BD into a giant envelope are those of Livio (1982), Livio & Soker (1983, 1984) and 
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Soker et al. (1984). For the majority of their calculations, the authors investigate whether 
the secondary gains or loses mass during the spiral-in, and how this affects the outcome of 
the CE process. Such an investigation is beyond the scope of the present paper. In those 
calculations where the authors assumed the secondary mass did not change, they concluded 
that the spiral-in would lead to merger in most cases. However, the calculations by Livio 
and Soker were one-dimensional, and thus crude compared with current three-dimensional 
hydro dynamical calculations (cf. Iben & Livio 1993; Taam & Sandquist 2000). 

In the following section, we describe the method used to model the present-day for- 
mation of CVs with BD secondaries. In section 3, we present results for the orbital period 
distribution of CVs forming with BD secondaries at the present time. We also describe 
the results from our investigation of the dependence of the formation rate of CVs with BD 
secondaries on the assumed value oi acE- In section 4, we present a discussion of the results 
and a comparison with observations. We close with a summary of our main conclusions. 

2. Method 

The population synthesis code developed by Politano (1988, 1996) to model the forma- 
tion of CVs has been modified for the purposes of this study. For a complete description of 
the code, we refer the reader to the original paper (Politano 1996), and here only discuss the 
modifications made. 

(1) The code developed by Politano (1988, 1996) employs a Jacobian technique to calcu- 
late the present-day CV birthrate and to map the ZAMS binary progenitor distributions to 
the ZACV distributions. The Jacobian approach provides a completely analytic framework 
from which to calculate theoretical statistics of close binary systems from the parent popu- 
lation of ZAMS binaries. However, in practice, it is difficult to test the effects of different 
values of certain input parameters (such as acE or the treatment of magnetic braking) on 
the computed ZACV population using the Jacobian code. A Monte Carlo technique was 
used by dcKool (1992) to compute the CV birthrate. He used identical or slightly modified 
physics as in Politano (1988, 1996) and found substantially the same results. A distinct 
advantage of the Monte Carlo approach over the Jacobian approach is that the Monte Carlo 
approach allows input parameters to be easily changed, and therefore lends itself nicely to 
calculating model populations for a wide range of input parameters. Since deKool (1992), 
several other authors have used a Monte Carlo technique in population synthesis calculations 
of CVs (e.g., di Stefano & Rappaport 1994; Han et al. 1995; Howell et al. 2001). We have 
therefore developed a Monte Carlo version of the Jacobian code, and have used the Monte 
Carlo code in the present study. The Monte Carlo code contains the same physics as the 
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Jacobian code, except as noted below, and it has been calibrated against the Jacobian code. 

(2) The birthrate calculations of Pohtano (1996) were modified to include a more realistic 
treatment of the radius of the BD as a function of mass and age. This treatment comes from 
coohng curves of BDs produced by I. Baraffe and co-workers (Baraffe et al. 1998, Baraffe 
2003, priv. comm.). We show some representative curves in Figure 1. The models in Figure 
1 and similar models (e.g., Burrows et al. 1997) indicate that BDs, as they cool, approach a 
radius that is ~ 0.1 R©. However, the cooling timescale to reach this radius can be as long 
as 10^ years for the most massive BDs, and the radius can decrease by a factor of four. Thus, 
unlike the calculations of Politano (1996), which assumed that the contraction phase to the 
main sequence for stars was short enough to be ignored, the cooling phase in BDs cannot 
be similarly ignored. To avoid any ambiguity, we note that the term "brown dwarf is used 
in this paper to denote an object whose mass is between 0.013 and 0.075 Mq (cf. Chabrier 
& Baraffe 2000). The lowest mass secondary considered in this paper is in accord with this 
definition. 

(3) In many population synthesis studies of CVs, the parent ZAMS binaries are assumed 
to be distributed over primary mass (Mp), mass ratio (q, defined as the ratio of the secondary 
mass to the primary mass), and orbital period (P) according to the following prescription: 

N{Mp, q, P) dMp dq dP = f{Mp) g{q) h(P) dMp dq dP (2) 

where f(Mp), g(q) and h(P) are distributions over the respective orbital parameters (e.g., 
Politano 1988, 1996; deKool 1992; di Stefano & Rappaport 1994; Han et al. 1995). In the 
present Monte Carlo code, the distribution of primary masses, f(Mp), is chosen from a Miller 
& Scalo (1979) initial mass function (IMF). The distribution of progenitor orbital periods, 
h(P), is assumed to be fiat in log P (cf. Abt 1983). Both of these choices are the same as 
those in Politano (1988, 1996). 

For a given primary mass, the secondary mass is chosen from the distribution, i(Ms) = 
Mpg(q). Therefore, the distribution of secondary masses is strongly infiuenced by the choice 
of g(q), the initial mass ratio distribution (IMRD) for the ZAMS progenitor binaries. The 
IMRD chosen by Politano (1988, 1996) was g(q) oc q. This gives a distribution of initial 
mass ratios that is strongly peaked toward q = 1 (i.e., equal primary and secondary masses). 
Observations of ZAMS spectroscopic binaries done after Politano's original work (Politano 
1988) suggest that the IMRD is not very strongly peaked toward q = 1, but is fiatter (i.e., 
g(q) oc q^, where a is very small; Duquennoy & Mayor 1991; Goldberg et al. 2003). For the 
calculation of the present-day orbital period distribution of ZACVs with BD secondaries, 
we have adopted an IMRD, g(q) dq = 1 dq. We note though that both Duquennoy & 
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Mayor (1991) and Goldberg et al. (2003) point out that there are large uncertainties in their 
derivation of the IMRD for very small mass ratios (q < 0.20). These uncertainties are due to 
the difficulty in observing the lines of the secondary star, since the light is dominated by the 
primary. Thus, we emphasize that the IMRD for extreme mass ratio binaries remains very 
poorly constrained and our choice, g(q) = 1, is by no means a definitive one. We further 
emphasize that, as in Politano (1988, 1996), we have assumed that the IMRD for binaries 
containing a ZAMS star and a BD is the same as the IMRD for binaries containing two 
ZAMS stars. This is an important assumption whose validity is discussed later in section 
4.2.3. 

In our Monte Carlo calculations, we begin with 10^ ZAMS binaries. This is suflicicnt to 
give a statistically significant sample of present-day ZACVs (~10^ systems). The age of the 
Galaxy is assumed to be 10^° yrs. To reproduce the Miller-Scalo IMF, we used Eggleton's 
(2001) Monte Carlo representation of it: 

Mp{u) = 0.19m[(1 - uf/^ + 0.032(1 - uy/*]-^ (3) 

where u is a random number uniformly distributed between and 1. As we are only 
concerned with the relative fraction of the total ZACV population that forms with BD 
secondaries, normalization of the above distributions has been ignored. 

3. Results 

The results of the population synthesis calculations are shown in Figures 2 and 3. In 
Figure 2, the present-day orbital period distribution of ZACVs is shown for the assumed 
choices, acE — 1 and g(q) = 1. The subset of ZACVs that contain BD secondaries is shown 
in the dashed histogram. Our model predicts that CVs that are presently forming with BD 
secondaries occupy an orbital period range from 46 min to 2.5 hrs. We further predict that 
18% of CVs at the present time arc forming with BD secondaries. Of these, the majority 
{'^ 80%) form below the orbital period minimum (78 min), with the remainder forming a 
tail reaching within the stellar domain. Thus, our model predicts that a significant fraction, 
~15%, of CVs are forming at orbital periods shorter than 78 minutes at the present time. 

In Figure 3, the fraction of present-day ZACVs that contain BD secondaries is plotted 
as a function of acE, the common envelope efficiency parameter, for 3 different IMRDs. The 
IMRDs chosen arc: g(q) oc q, g(q) oc 1, and g(q) oc q^"'^. Since the total number of ZACVs 
is sensitive to the choice of IMRD (cf. deKool 1992; Politano 1994), it may be true that the 
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fraction of CVs that form with BD secondaries also is sensitive to this choice. In hght of the 
aforementioned lack of strong observational constraints on the IMRD, we chose functional 
forms for g(q) that range from a distribution that strongly favors the formation of extreme 
mass ratio ZAMS binaries (g(q) oc q""^'^), to one that strongly favors the formation of ZAMS 
binaries that contain equal mass components (g(q) oc q). 

Not surprisingly, we find that the fraction of ZACVs with BD secondaries increases as 
the IMRD favors more and more extreme mass ratio progenitor binaries. For acE = 1, 
the value used in the calculation of the orbital period distribution in Figure 2, this fraction 
ranges from 8%, for g(q) cx q, to slightly greater than one-third (34%), for g(q) oc q~°'*^. For 
a given IMRD, the fraction of ZACVs with BD secondaries is relatively constant from acE = 
1.0 to acE = 0.3, and then begins to decrease for acE < 0.3. However, we find that even for 
very small values of the common envelope efficiency parameter {acE — 0.1), the fraction of 
CVs that form with BD secondaries is not zero. Rather, for acE — 0.1, this fraction ranges 
from ~ 3% to ~ 18%, depending on the IMRD chosen. 

4. Discussion 

4.1. Comparison with Observations 

Our models predict that the majority of CVs that form with BD secondaries populate 
an orbital period interval, 46 min < Porb < 78 min, that is not populated by CVs that formed 
with ZAMS secondary stars. This result, calculated using up-to-date BD models, confirms 
the previous result by Politano (1996) that ultra-short period CVs (Porb ^ 80 min) can be 
formed if the secondary is substellar at birth. Such systems therefore may provide a viable 
model for the observed ultrashort-period CVs, V485 Cen and RX 2329-1-06. We discuss this 
possibility below. In addition, there are a handful of CVs above the period minimum, in the 
orbital period range ~ 78-90 mins, that are strongly suspected of containing BD secondaries. 
Since, ZACVs with BD secondaries extend in orbital period to 2.5 hours, we also discuss 
below whether these handful of suspected BDCV systems may be CVs that formed with BD 
secondaries. 

We recognize that we are comparing a present-day orbital period distribution of ZACVs 
with BD secondaries with the observed, present-day orbital period distribution of CVs. 
Clearly, a more correct comparison would be between a secularly-evolved BDCV orbital 
period distribution and the observed one. However, the mass transfer rates in this particular 
sub-population of CVs are so low (~10~^^ Mq yr~^; Kolb & Baraffc 1999) that we do not 
expect the secularly-evolved distribution to be vastly different than the zero-age one. For 
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example, for the secular evolution of a 0.07 M© BD and a 0.6 Mq WD, Kolb & Baraffe (1999) 
found that it took 1.5 Gyrs for the system to evolve from 54 min to 78 min. Therefore, we 
do not anticipate that subsequent secular evolution of ZACVs with BD secondaries will 
evacuate the orbital period regime < 78 min. Nevertheless, our estimates of the fraction of 
ZACVs with BD secondaries should be regarded as upper limits. Plans to extend the present 
calculation to include subsequent secular evolution are in progress. 

4.1.1. V485 Cen and RX2329+06 

As mentioned earlier, V485 Cen and RX2329+06 are unusual CVs because their or- 
bital periods of 59 min and 64 min, respectively, are significantly shorter than the orbital 
period minimum for CVs. Suggestions as to why V485 Cen and RX2329+06 have peri- 
ods shorter than the period minimum have focused on the possibility that the secondary 
star in these systems had undergone significant nuclear evolution in its core prior to mass 
transfer (Augusteijn et al. 1996; Schenker & King 2002; Thorstensen et al. 2002; Uemura et 
al. 2002). Theoretical calculations by Schenker & King (2002) indicate that a CV contain- 
ing a secondary core in which 95% of the hydrogen has been burned into helium prior to 
mass transfer could reach a minimum period < 1 hr. A similar calculation by Thorstensen 
et al. (2002) found comparable results (see also Pylyser & Savonije 1988; Podsiadlowski, 
Rappaport & Pfahl 2002). 

While ZACVs with BD secondaries can account for the ultra-short orbital periods in 
V485 Cen and RX2329-I-06, recent spectroscopic observations of RX2329-I-06 by Thorstensen 
et al. (2002) have determined an unusually high temperature of ~ 4000 K for the secondary. 
Howell (2002, priv. comm.) has also found a similar temperature for the secondary in V485 
Cen. If these temperature determinations are correct, then this presents serious problems for 
our model. Even the youngest, hottest BDs do not have surface temperatures as high as 4000 
K (I. Baraffe 2002, priv. comm.; see also Chabrier & Baraffe 2000). It is possible that the 
close proximity of a hot white dwarf could cause significant heating of the BD atmosphere. 
Hauschildt (2002, priv. comm.) has performed preliminary calculations of the irradiation 
of a gas giant planet by a 15,000 K white dwarf at a distance of 1 AU. He did not find 
a significant increase in the temperature of the planet. However, the orbital separations in 
V485 Cen and RX2329-I-06 are much smaller, of order a few tenths of a solar radius. Harrison 
et al. (2003) have recently modeled the light curve of EF Eri, a CV with an orbital period 
of 81 min that is strongly suspected of containing a substellar secondary (see section 4.1.2. 
below) . They infer from their models that the secondary star in EF Eri is very cool (Tefr ~ 
900 K) and is irradiated by the WD. They estimate the irradiated side of the secondary has 
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a temperature, T^s ^ 1600 K. However, they note that proper modehng of the hght curve 
of EF Eri will require further work, such as the incorporation of BD atmospheres into the 
light curve program used by them. At the moment, therefore, it is not known whether the 
conditions that exist in V485 Ccn and RX 2329+06 could provide enough heating of the BD 
atmosphere to account for the observed temperatures. 

Other observed features in V485 Cen and RX2329+06 also favor the evolved secondary 
model. Both systems have Hel A6678/Hq; ratios that are higher than what is seen in typical 
SU UMa systems (Augusteijn et al. 1996; Thorstensen et al. 2002). This suggests that 
the accreted material is enriched in helium. In addition, the observational estimates of the 
secondary mass in V485 Ccn (M^ec > 0.14 M©; Augusteijn et al. 1996) and RX2329+06 
(Msec ~ 0.13 Mq, assuming a typical mass of 0.6 Mq for the WD; Thorstensen et al. 2002) 
are not consistent with a substellar secondary. 

4- 1-2. CVs above the Period Minimum Suspected to Contain BD Secondaries 

In a recent paper, Littlefair et al. (2003) reviewed the observational evidence for BD 
secondaries in CVs. They conclude that there is signficant indirect evidence for BD secondary 
stars in 6 CVs, all of which have orbital periods greater than 78 minutes. These systems 
are DI UMa (78.6 min), EF Eri (81.0 min), WZ Sge (81.6 min), AL Com (81.6 min), EG 
Cnc (86.3 min), and VY Aqr (90.8 min). As noted in the Introduction, according to the 
standard model of CV secular evolution (e.g.. King 1988), as the secondary star transfers 
mass to the WD, the orbital period of the system shrinks until, at a certain point, the 
secondary's thermal timescale becomes longer than the mass transfer timescale (Paczyhski 
& Sienkiewicz 1981; Rappaport, Joss & Webbink 1982). This point roughly coincides with 
the secondary having lost so much mass that it has become substellar (cf. Howell et al. 
1997). Further mass transfer then results in an increase in the orbital period of the system 
(Paczynski & Sienkiewicz 1981; Rappaport, Joss & Webbink 1982). From Figure 1, we see 
that the orbital period distribution of present-day ZACVs containing BD secondaries extends 
above 78 minutes to 2.5 hours. Therefore, at the present time, we cannot unambiguously 
say whether the above 6 systems are CVs that contain secondaries that were born substellar 
or became substellar as a result of mass loss during secular evolution. Below we compare the 
likelihood of each. 

Theoretical studies of the intrinsic (i.e., actual, not observed), current CV population 
predict that, according to standard CV secular evolution, 99% of the current population of 
CVs have orbital periods less than 3 hrs, and of these, 60-70% contain secondaries that have 
lost enough mass to become substellar (Kolb 1993; Howell et al. 1997, 2001). Under the 
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assumption of a constant CV birthrate throughout the Galaxy's history (e.g., Kolb 1993), 
our calculations predict an upper limit to the percentage of CVs in the period range 78 min-3 
hr, whose secondaries were born substellar, of 18%. This assumes all CVs that form below 
78 minutes quickly evolve to orbital periods greater than 78 minutes (as discussed above, 
this is clearly an unrealistic assumption but we use it as an extreme case to calculate an 
upper limit). Therefore, assuming standard secular evolution, it is at least ~ 4 times more 
likely that CVs in the period range 78 min-3 hr contain secondaries that became substellar 
due to mass loss rather than were born substellar. 

We note that recent observations of the rotation rates of low-mass stars in open clusters 
(Terndrup et al. 2000) have raised some concerns about the validity of the standard model 
of CV secular evolution. An empirical angular momentum loss prescription inferred from 
these observations (Andronov et al. 2003) suggests that the strength of the mean magnetic 
braking in CVs may be on the order of 100 times less than that assumed in the standard 
model for CVs above the period gap (Porb > 3 hrs; cf. King 1988). Such a drastic reduction 
in the mean magnetic braking would imply proportionately longer lifetimes for CVs above 
the gap. If the magnetic braking above the gap is indeed as low as suggested by Andronov 
et al. (2003), then King & Schenker (2002) have postulated that the 78 min cutoff in the 
orbital period distribution of CVs may be a feature dictated by the age of the Galaxy. 
That is, 78 minutes may simply be the shortest orbital period achievable by a CV within 
the age of the Galaxy, given the much longer evolutionary timescales above the gap due 
to reduced magnetic braking. If true, this would drastically reduce the large population of 
CVs (predicted by the standard model) whose secondaries became substellar as a result of 
mass loss, and thereby increase the likehhood that systems like WZ Sge contained substellar 
secondaries at birth. 

The rotation rate data from open clusters also indicate no significant reduction in mag- 
netic braking once a star becomes fully convective (Terndrup et al. 2000) , as assumed in the 
standard model (e.g., Rappaport, Verbunt & Joss 1983). Thus, magnetic braking should 
also operate below the period gap, whereas in the standard model, the only angular momen- 
tum mechanism driving mass transfer below the gap is gravitational radiation (Paczyhski 
& Sienkiewicz 1981; Rappaport, Joss & Webbink 1982). Andronov et al. (2003) infer an 
angular momentum loss rate below the gap that is ~ 1.5 times higher than gravitational 
radiation alone. Such higher mass transfer rates below the gap could completely exhaust 
the mass of the secondary within the age of the Galaxy, and thus significantly reduce the 
number of secondaries that became substellar as a result of mass loss (cf. Patterson 1998). 
It should be noted that the empirical angular momentum loss prescription by Andronov et 
al. (2003) does not provide a satisfactory explanation for the existence of the period gap in 
CVs, or the observed brightnesses of CVs above the gap. 
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4.2. Where are All the CVs with Orbital Periods less than 78 mins? 

There are approximately 500 CVs with known or suspected orbital periods in the Ritter 
& Kolb (2003) catalog. Regardless of whether V485 Cen or RX2329+06 contained an evolved 
secondary or a substellar secondary at birth, the observed number of CVs in the period range 
46 min-78 min is still considerably less than predicted by our models of ZACVs with BD 
secondaries (15% or ~ 75 CVs) for an IMRD, g(q) = 1, and acE — 1- Even with the 
recognition that 15% is an upper limit, the discrepancy between our predicted number of 
CVs with orbital periods less than 78 mins and the observed number is significant. In the 
sections below, we discuss possible reasons for the cause of this discrepancy. 



4-2.1. Observational Selection Effects 

The mass transfer rates in ZACVs with BD secondaries are predicted to be low (~ 10~^^ 
M0 yr~^; Kolb Sz Baraffe 1999), therefore these systems should be faint (m„ > 20; Politano 
et al. 1998) and difficult to observe. Consequently, it is possible that these systems exist, but 
are too faint to be observed currently in quiescence. Indeed, the 6 systems suspected of con- 
taining BD secondaries discussed above were discovered as a result of their outbursts. The 
recurrence times betweeen outbursts in these six systems is longer than most dwarf novae 
(cf. Warner 1995), and this may be a property of CVs that contain substellar secondaries. 
Presumably, the longer recurrence times are a result of the comparatively lower mass transfer 
rates in these systems than typical dwarf novae. Kolb & Baraffe (1999) calculated evolu- 
tionary tracks for CVs with initial secondary masses between 0.04 M0 and 0.21 Mq. They 
computed mass transfer rate as a function of orbital period, and found that the evolutionary 
tracks merge into a common track for orbital periods above approximately 70 minutes, re- 
gardless of whether the secondary star was initially substellar or became substellar as a result 
of mass loss (cf. Fig. 1 in Kolb & Baraffe 1999; we note that the BD models used in Kolb 
& Baraffe (1999) are very similar to the ones used in the current paper.) Their tracks also 
show that a CV that formed with a substellar secondary below 70 mins (0.05 < Msec < 
0.075 Mq) has a mass transfer rate higher than its evolved counterpart above the period 
minimum whose secondary became substellar as a result of mass loss (again, cf. Figure 1 in 
Kolb & Baraffe 1999). Thus, although observational selection effects cannot be ruled out, 
if we are observing outbursts in systems with presumably substellar secondaries above the 
period minimum (e.g., WZ Sge, AL Com, etc.), then it is difficult to understand why we 
would not be able to observe outbursts in systems formed below the period minimum, if the 
mass transfer rates are the same or higher. 
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4-2.2. Common Envelope Efficiency Parameter for Very Low Mass Secondaries 

Our results indicate that CVs will form with BD secondaries unless acE is extremely 
small (< 0.1). At first glance, this is a surprising result. However, in the population synthesis 
code used in these calculations, the orbital separation at the end of the CE phase is only 
limited by the condition that the secondary does not fill its Roche lobe during the spiral-in 
process. This effectively allows spiral-in for such small objects to an orbital separation of 
~ a few tenths of a solar radius. Thus, the decrease in the orbital energy reservoir due to 
such small secondary masses is offset by the very small orbital separations these systems 
may achieve upon exit of the CE phase. Therefore, only if the transfer of orbital energy into 
the envelope is extremely inefficient will merger occur. 

Such extremely small values of acE for the entire population of CV progenitors would 
produce a total number of CVs in the Galaxy that would be difficult to reconcile with even 
the lowest estimates of the CV space density (~10^^pc^'^; Patterson 1984). In addition, 
estimates of the efficiency parameter extrapolated from detailed hydrodynamical calculations 
of the CE phase have typically indicated acE ^ 0.3 (e.g., Iben & Livio 1993; Taam & 
Sandquist 2000). 

It may be that acE is a function of the secondary mass. In particular, if acE decreases 
sharply as the secondary mass approaches the substellar regime, this could account for the 
lack of observed CVs below 78 minutes. For example, it may be true that the luminosity 
generated by the drag torque on the secondary becomes comparable to the luminosity of the 
giant star for very low mass secondaries. In this case, merger cannot be avoided since the 
time scale for the generation of frictional energy would be too long to dynamically expel the 
envelope (R. Taam 2002, priv. comm.; see also Iben & Livio 1993 and references therein). 

A second possibihty is that acE is different for CE evolution involving primaries that 
fill their Roche lobes on their initial ascent of the giant branch (i.e., before He ignition) than 
for primaries that fill their Roche lobes as AGB stars. There have been several detailed 
hydrodynamical calculations investigating the spiral-in of a low mass secondary star into an 
AGB primary (cf. Iben & Livio 1993). However, there has been only one detailed study of the 
spiral-in of a low-mass secondary into an RGB primary (Sandquist et al. 2000). It appears 
from this study that acE could be different for CE evolution involving an RGB primary than 
for an AGB primary. However, further detailed study of this scenario is necessary before 
any firm claim can be made. Population synthesis calculations of the formation of CVs that 
include a variable acE are currently underway (Politano & Quandt 2003). 

Finally, although our results indicate that there is a sufficient energy reservoir in the 
orbit to unbind the envelope even if acE is quite small, Taam and co-workers have noted 
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that the drag torque on the secondary must also approach zero, or else the secondary will 
continue to spiral in towards the core. They find that this dynamical constraint on the 
drag torque requires the presence of steep density gradients within the common envelope 
(Yorke, Bodenheimer & Taam 1995; Terman & Taam 1996). To the best of our knowledge, 
no population synthesis code for CVs (including our own) includes this dynamical constraint 
on the drag torque in its treatment of the common envelope phase. 

4-2.3. IMF and Mass Ratio Distribution for Binaries with Suhstellar Companions 

Several studies have shown that the daughter ZACV population is strongly dependent 
on the assumed IMRD chosen for the parent ZAMS binary population (e.g., deKool 1992; 
Politano 1994). As evident in Fig. 3, the population of ZACVs with BD secondaries is 
similarly strongly affected by the assumed choice of the progenitor IMRD. In the present 
study, we have assumed that the IMRD for binaries containing a ZAMS star and a BD is the 
same as the IMRD for binaries containing two ZAMS stars. As described below, preliminary 
results for the BD mass function and the frequency of BDs in binaries suggest that this 
assumption likely is not valid. 

Observational searches for BDs in young clusters (Martin et al. 1998, 2000; Bouvier et 
al. 1998; Hambly et al. 1999), in star forming clusters (e.g., Luhman 2000; Luhman et al. 
2000; Bouvier et al. 2001), and in the field (Reid et al. 1999) suggest that the mass function 
for single BDs is approximately flat or rises slowly from the substcllar regime to ~ 0.6 M©, 
and then rolls over into a power law from ~ 1 M© to higher masses. Estimates of the number 
density of nearby BDs from surveys such as the Deep Near-Infrared Survey (DENIS; Delfosse 
et al. 1999) and the Two Micron All-Sky Survey (2MASS; Kirkpatrick et al. 1999; Burgasser 
et al. 1999) are consistent with such a BD mass function, although clearly we are dealing 
with small number statistics. 

Investigation of the mass function of BDs in binaries is an even more recent undertaking. 
Observational data in the past few years have allowed some preliminary statements to be 
made. These data indicate that BD companions to F-M stars in "close" orbits (< 3 AU) are 
rare (the so-called "brown dwarf desert;" cf. Marcy & Butler 2000). Recent observations by 
Gizis et al. (2001) indicate that no such desert exists for BD companions to F-M stars in 
wide orbits (> 1000 AU). 

From the present calculations, we find that CVs that are currently forming with BD 
secondaries have evolved from a parent, ZAMS population of binaries with orbital separations 
< 3 AU and primary star masses in the range ~ 1-10 Mq. In Figure 4, we present histograms 
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of the number of progenitors as a function of (a) the progenitor ZAMS primary mass and (b) 
the progenitor ZAMS orbital separation, respectively. Inspection of Figure 4a indicates that 
~ 75% of the primaries in the parent, ZAMS population have masses less than ~ 1.6 Mq, 
and thus have spectral types F or later. Inspection of Figure 4b shows that the vast majority 
of the progenitors have orbital separations less than 1 AU. Further, (although this cannot 
be discerned in the figure), the largest orbital separation in the progenitor population is 2.98 
AU. Interestingly, these ranges in orbital separation and primary mass for approximately 
three-fourths of the progenitors of ZACVs with BD secondaries lie directly within the brown 
dwarf desert. It is possible, therefore, that the lack of observed CVs with orbital periods 
less than 78 minutes may be due to a lack of progenitor binaries from which to draw. 
Although the existence of the brown dwarf desert seems fairly well established, we believe 
that further detailed observations of substellar objects in close binaries are needed before 
any firm conclusions can be drawn. Furthermore, our understanding of acE for CE evolution 
involving very low-mass secondaries needs to improve before we can separate its effect on 
the formation of ZACVs with BD secondaries from the effect of the IMRD. 

The result that all progenitors of ZACVs with BD secondaries have ZAMS orbital sep- 
arations less than 3 AU is nevertheless quite significant, and underlines the potential im- 
portance of observations of CVs with ultra-short orbital periods. Assuming the effect of 
acE on the results can be isolated, then it may be possible to use observations of CVs with 
orbital periods less than 78 mins, along with the models presented in this paper, to provide 
useful information about the progenitor population of ZAMS -|- BD binaries. For exam- 
ple, by comparing our models with these observations, we may be able to place constraints 
on the BD mass function and/or IMRD. These constraints would nicely complement, and 
possibly augment, our understanding of the BD mass function/IMRD in MS binaries from 
observations of these systems in clusters and in the field, for example. In addition, a firm 
null result for the detection of CVs with BD secondaries and orbital periods less than 78 
min has the potential to provide strong, independent evidence for the existence of the brown 
dwarf desert. 

We close this section with a few speculative remarks regarding the formation of BDs in 
binaries. According to standard stellar evolution theory, primary stars with masses < 1 Mq 
will not evolve off of the main sequence within the age of the Galaxy. One might therefore 
expect to see a relative increase in the number of BD-MS star binaries with primary masses 
below 1 Mq. However, the brown dwarf desert extends in primary mass well below 1 Mq to 
M-type primaries. Secondly, an observational point that has not received as much attention 
as the brown dwarf desert is that there is no observed scarcity of L dwarf-L dwarf binaries 
(cf. Reid et al. 2001; Chabrier 2002 and references therein). These statements suggest that 
there is something different about the frequency or formation rate of BD-MS star binaries 



-16- 



compared with MS-MS binaries or BD-BD binaries. Much attention has been paid to the 
BD mass function in binaries. We suggest that perhaps the more important distribution is 
the mass ratio distribution in binaries containing BDs. The existing observational data are 
also consistent with the following statement: BDs in binaries with approximately equal mass 
ratios are not rare, while BDs in binaries with extreme mass ratios arc, at least at separations 
< 3 AU. This raises some interesting questions: Is the frequency of BDs in binaries strongly 
affected by the mass ratio? Is the mass ratio distribution in extreme mass ratio binaries 
independent of the total mass of the binary? For example, is there a corresponding relative 
scarcity of ZAMS binaries with extreme mass ratios compared to those with nearly equal 
mass ratios? Goldberg et al. (2003) have recently derived the mass ratio distribution for a 
sample of 129 spectroscopic binaries with orbital periods less than 2500 days. They find that 
there is a sharp drop in the mass ratio distribution for mass ratios less than ~ 0.2 (again, 
secondary mass/primary mass). However, they caution that there are large uncertainties in 
their derivation for such small ratios, and they are hesitant to ascribe a statistical significance 
to this finding. We believe each of these questions merit further study. 



5. Summary and Conclusions 

We summarize our key findings and conclusions below: 

(1) CVs that are presently forming with BD secondaries are predicted to occupy an 
orbital period range of 46 min to 2.5 hrs. This result is in substantial agreement with the 
previous result by Politano (1996). For the assumed choices of acE — 1 and a mass ratio 
distribution in the ZAMS progenitor binaries, g(q) dq = 1 dq, we predict that 18% of all 
present-day ZACVs contain secondaries that are BDs. Of this 18%, the majority (~ 80%) 
form below the orbital period minimum (78 min), with the remainder forming a tail reaching 
within the stellar domain. Thus, our model predicts that a significant fraction, ~15%, of 
CVs are forming with orbital periods shorter than 78 minutes at the present time. 

(2) CVs that form with BD secondaries could provide a viable model to explain the 
ultra-short periods in the unusual CVs, V485 Cen (Porb = 59 min) and RX2329+06 (Porb = 
64 min). However, because of the unusually high temperatures derived for the secondaries 
(~ 4000 K) in these systems (Thorstensen et al. 2002; Howell 2002, priv. comm.), it appears 
unhkely that these systems contain a BD secondary, unless irradiation of the BD by the 
nearby WD has a significant effect on the surface temperature of the BD. These high tem- 
peratures for the secondary are more easily explained with a model in which the secondaries 
in V485 Cen and RX2329+06 are the cores of evolved MS stars (Schcnkcr & King 2002; 
Thorstensen et al. 2002; Uemura et al. 2002). This model is also in better agreement with 
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the estimated secondary masses and enhanced He/H ratio observed in these systems. 

(3) CVs that form with BD secondaries could also provide a model for six CVs above the 
period minimum, DI UMa, EF Eri, WZ Sge, AL Com, EG Cnc, and VY Aqr, for which there 
is signficant indirect evidence that the secondary stars in these systems are BDs (Littlefair et 
al. 2003). These six systems have orbital periods between ~ 78-90 min. However, according 
to the standard model of CV secular evolution, it is at least four times more likely that the 
secondaries in these six systems became substellar as a result of mass loss, rather than were 
born substellar. It is important to note though that the validity of the magnetic braking 
prescription used in the standard model has been recently called into question by observations 
of the rotation rates of low-mass stars in open clusters (Terndrup et al. 2000; Pinsonncault, 
Andronov & Sills 2002). A prescription for magnetic braking inferred from the open cluster 
data (Andronov et al. 2003) could imply a significant reduction in the number of CVs whose 
secondaries became substellar as a result of mass loss (cf. Patterson 1998; King & Schenker 
2002). If true, this would correspondingly increase the hkehhood that WZ Sge, AL Com, 
etc. are CVs that formed with BD secondaries. 

(4) Assuming a constant CV birthrate throughout the Galaxy's history, the fraction of 
CVs that are presently forming below the period minimum according to our model, 15%, 
represents an upper limit to the predicted fraction of all CVs with orbital periods less than 78 
minutes at the present time. Of the ~ 500 CVs with known orbital periods, 15% represents 
75 systems. The number of CVs with orbital periods less than 78 minutes is currently two, 
V485 Cen and RX 2329+06, and it is dubious whether these systems contained substellar 
secondaries at birth. Thus, even with the recognition that 75 systems is an upper limit, the 
observed number of CVs in the period range 46 min-78 min is still considerably less than 
the number predicted by our models. It does not seem likely that the discrepancy is due 
to observational selection effects. The mass transfer rates and, hence, recurrence times in 
CVs with BD secondaries below the period minimum should be very similar to systems like 
WZ Sge and AL Com (cf. Kolb & Baraffe 1999). Consequently, it is difficult to understand 
why we would not have observed outbursts in them. Ongoing searches for faint CVs in 
observational surveys such as 2MASS (Hoard et al. 2002; Wachter et al. 2003), the Sloan 
Digital Sky Survey (SDSS; Szkody et al. 2002) and the Faint Sky Variability Survey (FSVS; 
Howell et al. 2002) may be able to eventually detect such systems in quiescence, assuming 
they exist. 

(5) We have examined the dependence of the present-day formation of CVs with BD 
secondaries on the assumed value of the common envelope efficiency parameter, acE- We 
find that even for values of acE as low as 0.1, CVs with BD secondaries are still able to be 
formed. For acE — 0.1, our calculations predict that the present-day fraction of ZACVs with 
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BD secondaries is 8% for an fiat IMRD, g(q) dq = 1 dq. For an IMRD that strongly favors 
tlic formation of binaries with mass ratios near unity, g(q) oc q, this fraction is decreased to 
3%. For an IMRD that strongly favors the formation of binaries with extreme mass ratios 
near, g(q) oc q~°'^, this fraction is increased to 18%. A firm null result for the detection 
of CVs with BD secondaries having orbital periods less than 78 mins could therefore imply 
that the transfer of orbital energy into the common envelope is extremely inefficient for 
systems with very low-mass secondaries. Further, as it is unlikely that the CE process is 
that inefficient for the entire CV population, such a null detection would also then imply 
that the efficiency of the CE process is a function of secondary mass. If true, it would then 
seem extremely coincidental if acE decreases abruptly exactly at the substellar transition 
mass. Rather, one would expect that the decrease is smoother and begins in the stellar 
domain. Consequently, determining if there is a cutoff mass for the secondary, below which 
merger is inescapable, may have important implications for our theoretical understanding of 
the orbital period distribution in CVs below, and possibly within, the period gap. 

(6) Lastly, CVs that form with BD secondaries have evolved from progenitor ZAMS + 
BD binaries with primary masses in the range ~ 1-10 M0 and with orbital separations < 
3 AU, for an assumed IMRD that is flat (g(q) = 1). Furthermore, 75% of these progenitor 
systems have primary masses < 1.6 Mq, and therefore spectral types F or later. The orbital 
parameter ranges in primary mass and orbital separation for these 75% coincide with the 
so-called "brown dwarf desert." Comparison of the models presented in this paper with ob- 
servations of ultra-short period CVs may therefore be able to provide important information 
regarding the formation and distribution of BD -|- MS binaries with orbital separations less 
than 3 AU. 

In summary, we believe that the study of the formation of CVs with BD secondaries 
can provide useful insights concerning our understanding of common envelope evolution in 
close binary systems with very low mass secondaries, and the distribution of BDs in MS 
binaries having orbital separations < 3 AU. We urge theorists who are performing detailed 
hydrodynamical calculations of the CE phase to investigate in more detail the spiral-in of very 
low-mass objects in order to determine if there is indeed a cutoff mass below which merger 
is unavoidable. We are extremely grateful to the increasing numbers of observers involved 
in searches for faint CVs. These searches ultimately hold the key to our understanding 
of the global CV population. We urge observers involved in such searches, especially in 
the infrared, to be particularly vigilant regarding ultrashort-period CVs (~ 40-80 mins). 
For the reasons discussed in this paper, these systems may provide important guidance for 
constructing models of CV formation and secular evolution. This is particularly important at 
this time, as we reexamine, and in some cases abandon, long-held beliefs about the formation 
and secular evolution of CVs. 
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Fig. 1. — Radius as a function of time for BDs and low mass stars of solar composition. 
Numbers labeling curves indicate mass in solar masses. The time is measured in years from 
the formation of the BD. These coohng curves were taken from detailed model calculations 
by Baraffe and co-workers (Baraffe et al. 1998; Baraffe 2003, priv. comm.) 
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Fig. 2. — The predicted distribution of orbital periods in present-day ZACVs for the assumed 
choices: acE = and g(q) dq = 1 dq. The distribution has been divided into two histograms. 
The dashed histogram shows systems that contain BD secondaries. The sohd histogram 
shows systems that contain stellar secondaries. The total distribution is the sum of these 
two histograms. 
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Fig. 3. — The fraction of prcscnt-day ZACVs that contain a BD secondary as a function of 
acE for three different IMRDs. The dashed hne is for g(q) oc q~'^'^, the sohd hne for g(q) oc 
1, and the dashed-dotted hne for g(q) oc q. 



-27- 




02463 10 1 2 3 

ZAMS progenitor primary mass (M^) ZAMS progenitor orbital separation (AU) 



Fig. 4. — (a) The distribution of ZAMS primary masses in the progenitors of CVs that are 
forming with BD secondaries at the present epoch, (b) The distribution of ZAMS orbital 
separations in the progenitors of CVs that are forming with BD secondaries at the present 
epoch. The binaries in these distributions arc the parent population of the ZACV binaries 
shown in the dashed histogram (BD secondaries) in Figure 2. (Wc note though that we 
increased the number of initial binaries in our Monte Carlo calculations from 10^ to 10^ to 
ensure reasonable statistics for the entire range of primary masses in the above distribution.) 



